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Properties of scalar mesons below 1.0 GeV have been studied by regarding them as 
hadronic molecular states. Using the effective Lagrangian approach, we have calculated their 
leptonic decays, strong decays and productions via the (p meson radiative decays. Comparing 
our results with that given in the literature and the data, we conclude that it is difficult to 
arrange them in the same nonet if some of them are pure hadronic molecules. 



I. INTRODUCTION 



Understanding the nature of scalar mesons is a prominent topic in the past 30-40 years. The 
importance of the nature of the scalar mesons is that, because the properties of the scalar mesons, 
especially that with masses below 1.0 GeV, are difficult to be understood in the constituent quark 
models, the study of the light scalar mesons can help us to understand the nonperturbative proper- 
ties of QCD. Moreover, because the scalar mesons have the same quantum numbers as the vacuum, 
it can help us to reveal the mechanism of symmetry breaking which is, up to now, one of the most 
profound problems in particle physics. 

Many properties of scalar mesons are not so clear although they have been investigated for 
several decades. Experimentally, it is difficult to identify the scalar mesons because of their large 
decay widths which cause strong overlaps between resonances and grounds and also because several 
decay channels open up within a short mass interval. And due to these problems many data on 
scalar mesons are not so precise, so that it is not so easy to reveal their underlying structures [l|,y,0]. 
Theoretically, there are many calculations based on different models, but it seems that one cannot 
rule out some of them based on the present measurements. Currently, the observation shows that 
the known 0^^ mesons below 2.0 GeV can be classified into to classes: One class with masses 
below (or near) 1.0 GeV and the other class with masses above 1.0 GeV [l|. 

To study light scalar mesons, one problem we have to confront is how to classify the present 
observed scalar objects. One opinion is that the scalar objects with masses below 1.0 GeV, including 
two isosinglets (t(600) arid ^(980), one isotriplet ao(980) and two isodoublets KQ{SdQ), can be 

7, ^, y, [3]. On the contrary, another opinion, inspired by linear 



classified into one nonet 



2 



sigma model and unitary quark model, is that (t(600), /o(980), ao(980) and -fCg (1430) form a scalar 
nonet Q, 3]. 

Concerning the structures of light scalar mesons, they are still open questions so far although 
many attempts have been made to understand them in the literature. For the light scalars mesons 
with masses below 1.0 GeV, some people, by considering their dominant two-body decays, believe 



that they are multiquark (or multiquark dominant) states 



molecular states (ij. 



17, 



18|, 
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l4, ll5|,Q or hadronic 



211 ] ( Note that some references list here do not different 



multiquak state from molecular state.). Alternative^ 



mesons were also investigated in the qq picture 
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aroperties of some of these light scalar 



24( 1 . Moreover, in Ref. 25|, the spectrum 



of light scalar mesons below 1.0 GeV were studied in the qq picture with including the instanton 
effect. And it was found that the qq components with the instanton or tetraquark effect can explain 
the spectrum at a certain level. In Ref. {2^, we also studied the the effect of instanton-induced 
interaction in light meson spectrum on the basis of the phenomenological harmonic models for 
quarks. For the light scalar mesons with masses above 1.0 GeV, the potential model calculation 
indicates that they are qq states. However, concerning the spectrum of the scalar mesons with 
masses above 1.0 GeV, since there are three isoscalar states, /i(1370), /o(1500) and /o(1710), one 
believes that there is a glueball candidate among them. Based on the recent lattice results, the 
mixing between glueball and qq components in these three objects were fit in Ref. j^]- It should be 
noted that, the two problems we mentioned above are not independent. If we know the classification 
of the scalar mesons and the structures of some states, we may deduce the structures of other states 
according to the classification, and inversely, if the structures of the scalar mesons are confirmed 
the classification becomes obvious. 

In this paper, since some scalar mesons (mainly /o(980) and ao(980)) have been studied in the 
literature and the numerical results yielded there are consistent with the data, we will study the 
properties of all the scalar mesons with masses below 1.0 GeV by regarding all of them as pure 
hadronic molecules and check whether they form one nonet in the effective Lagrangian approach. 
Our logic is, if the scalar meson with masses below 1.0 GeV form a nonet, they should have the 



same structures. Anc. 
literature 



17 



181, 



3, 



if some of them can be interpreted as hadronic molecules as discussed in the 



20 



21[, all the elements in the same nonet should be interpreted as hadronic 



molecules. If our start point is reasonable, the yielded results should be consistent with the data. 
Otherwise, it is difficult to classify them into the same nonet, at least in the the hadronic molecular 
interpretation. 

In the molecular picture the interaction of scalar meson to its constituents can be described by 
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the effective Lagrangian. The corresponding effective coupling constant is determined by the 
compositeness condition Z = which was earlier used by nuclear physicists and is being widely used 



by particle physicists (see the references in 



28|] V In Refs. 



2S 



to study the newly observed charmed mesons 



28 



28, 



301 1 this method has been applied 



30| and the decay properties calculated there 



are consistent with the observed data. We had employed the above technique to predict the decay 
properties of the bottom-stran ge m esonst31l] and recently we applied this method to studied the 



baryonium picture of X(1835) 



have been studied in Refs. 



13, 



32 



18 



. The production and decay properties of some scalar mesons 



19,^ 



21[ with the help of this technique by regarding them as 



pure hadronic molecules and the numerical results yielded there agree with the data. 

It should be mentioned that, in the tetraquark interpretation of the scalar mesons with masses 
below 1.0 GeV, the diquark (qq) is in the 3 representation {us,ds,ud) of SU{3) group so there 
are at most two valence strange quarks in the (j(600) and /o(980) as shown in Fig. [TJ However 
in the hadronic molecular interpretation, an isosinglet r] constituent in addition to the isodoublet 
(K~^,K^) should be introduced to form a nonet. So that there is r]r] component, or {ss){ss) at 
the quark level, in o"(600) and /o(980) which makes it is difficult to understand the scalar meson 
spectrum from the quark level. Fortunately, in our present work we adopt the scalar meson 
spectrum as input, with the help of compositeness condition, to determine the coupling constant 
Qg between the molecule and its constituents. 

In our explicit calculation of the effective coupling constant g^, we will consider two cases, 
the one is local interaction case and the other one is nonlocal interaction case with including a 
correlation function which describes the distribution of the constituents among the molecules. As 
in our previous works, we will introduce a typical scale parameter to describe the finite size of 
the molecules. The numerical results show, in the region where the coupling constant g^ is stable 
against As, the yielded value of g^ for the corresponding molecular scalar meson, is consistent with 
that yielded from the local interaction picture. So that in the other calculations we will take the 
local interaction vertex. For other interactions, we will resort to the phenomenological Lagrangian 
and borrow the relevant coupling constants from the existing literature. 

To determine the mixing angle 9s between cj(600) and /o(980), we will apply the data for 
/o(980) 2j. It is natural that only with this data one cannot fix 9s uniquely and to fix 9s a well 
measured process of o"(600), for example (t(600) — > 27, is necessary. Concerning that the coupling 
constant g^ depends on the mass of cj(600) closely, we will not use this data, and will not discuss 
the physics of cr(600) in the present work. 

With the determined coupling constant g^ we will calculate leptonic decay constants of the scalar 
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mesons and the relevant leptonic decay widths. And from the yielded numerical results, we find 
that the leptonic decay widths of the scalar mesons are too small to be measured at present. With 
our hadronic interpretation, the strong decays of the scalar mesons to two-pseudoscalar mesons 
have also been calculated. The numerical results for the strong decays indicate that it is difficult 
to arrange the scalar mesons with masses below 1.0 GeV into the same nonet in the hadronic 
molecular interpretation. To study the productions of the scalar mesons /o(980) and ao(980) in 
the radiative decays of (j) meson, we will includ the final state interaction effect due to the two 
axial- vector mesons nonets with quantum numbers J^*^ = and 1^ . We find the final state 
interaction plays a negligible role in the production rates. 

This paper is organized as the following: In section [Til we discuss the molecular structures of the 
light scalar nonet with masses below 1.0 GeV and calculate the coupling constant g^, the leptonic 
decay constants and the leptonic decay widths. In section IIIH the strong decays of the scalar 
mesons are calculated in the framework of effective Lagrangian approach. The productions of the 
/o and oo mesons in the radiative decays of (j) meson are studied in section IIVI Our discussions 
and conclusions are given in the last section. 

II. MOLECULAR STRUCTURES OF THE SCALAR MESON NONET BELOW 1.0 GEV 

In this section, we will discuss the properties of the scalar meson nonet below 1.0 GeV in 
the hadronic molecular explanation, i.e., regarding them as two-pseudoscalar-meson bound states. 
At first, we will identify pseudoscalar meson contents of the scalar mesons and construct the 
effective Lagrangian describing the interaction between scalar meson and its constituents. Then 
the magnitude of the coupling constant is calculated with the help of compositeness condition 
and the leptonic decay constants and the leptonic decay widths of the scalar mesons are yielded 
by the standard loop integral. 

A. Molecular Structures of the Scalar Nonet 

Now, we are in the position to identify the constituents of light scalar meson nonet as hadronic 
molecules. As we pointed above that we will accept the opinion that the light scalar mesons: two 
isosinglets a (600) and /o(980), one isotriplet ao(980) and two isodoublets i^o(800), can be classified 
into one nonet, we should have three elements at least to form this nonet. And considering the 
physical requirement that the bound state should have a mass smaller than the threshold of the 
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corresponding constituents, the constituents of the hght scalar mesons should be K~^,K^ and r/. 
With these considerations, we may arrange the two-pseudoscalar-meson bound states in the (I3, Y) 
plane (where Y is the hypercharge which is relative to the strangeness S via relation Y = S for 
mesons and relates the charge Q and the third component of isospin via relation Q = l3 + (l/2)y) 
as Fig.[TJ As a comparison, we also illustrate the quark contents of the scalar nonet in the tetraquark 
interpretation in Fig. [H 
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FIG. 1: Hadronic molecules (left) and tetraquark explanations (right) of scalar mesons in the {I^^Y) plane. 

Comparing with the quantum numbers of the light scalar mesons, we write down the hadronic 
molecular structures of two isodoublets and isotriplet as 

|i^*+(800))\ ^ y ( l^o°(800)) 

|i^o*°(800)) j V^V' V 1^0" (800)) 

2ft (980)) \ ( K+K^\ „ 1 ^ _^ ^ 

° =1 ; |a0(980)) = -=(|K+K-)-|K0i^0)) (1) 

2o(980)) J \k^K- J V2 

and two singlets as 

58 = -^{\K+K-) + \K'K')-2\r^7^)) 

So = -^{\K+K-) + \K^K') + \rjr^)) (2) 

The physical states o"(600) and /o(980) should be mixing states of Ss and 5o via the mixing angle 
defined as 



/o(980) \ _ I cos 05 sin 6s 
cr(600) / I -sin 6*5 cos 6*5 









\ So / 



(3) 
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In terms of the constituents and the mixing angle 9s, the physical (t(600) and /o(980) mesons can 
be expressed as 



/o(980) 
a(600) 



(- 



sin 9s cos 9s 



,sm 9s 2cos( 



){\K+K-) + \K'>K'>)) + { 



V6 

, cos 9 s 2 sin 9 



In the case of ideal mixing with cos 9s = and sin 9s = \j\ one has 



+ 



(4) 



/o(980) 



{\K+K-) + \K^K^))] cj(600) = |r?r?) 



V2 



(5) 



This is the case in which the hadronic molecular picture of /(980) was studied in some literature, 



for example Refs. 



2C 



2l| . In our present work, we will take the mixing angle as a parameter 



and fit it from the data in the following. 

In summary, with the above discussions, one can write the scalar meson nonet in a matrix form 

as 

/ a[](980) + + ^i^So ^/2a+(980) ^27^*^(800) ^ 



1 

72 



V 



\/2ao (980) 



(980) + + -^Sq ^/2K*°(800) 



(6) 



And the ideal mixing reduces to 



'ideal 



1 

71 



( aO(980) + /o(980) 



\/2a(j"(980) \/2K(;+(800) 
](980) + /o(980) ^/2Ko*°(800) 



\ 



(7) 



^/2ao (980) 

\^ ^2^0^(800) ^Kl^{m^) ^/2ct(600) / 

It should be noted that in the hadronic molecular interpretation there are more quark contents 
in the relevant scalars compared to the tetraquark interpretation as shown in Fig. [TJ This makes 
it seems that it is difficult to understand the scalar meson spectrum in the molecular picture. 
However, at present, we do not attempt to calculate the scalar meson spectrum but take the scalar 
meson masses as input to determine the coupling constants. 

The effective Lagrangian, without including the distribution of the constituents among the 
molecules, describing the interaction between the scalar mesons and their constituents can be 
written as 



(8) 
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where the upper index "LC" denotes the local interaction vertex, S denotes the scalar meson matrix 
which was give in Eq. Q and P denotes the constituent pseudoscalar meson matrix 



P 



(9) 



V ^ J 

To include the distribution function which illustrates the distribution of the constituents in the 
molecules, one should modify the interaction Lagrangian ([8]) as 

^"s"- = J d^y^s{.y^)P}{x + u,y)S,,[x)P,{x - u^^y) (10) 

where the upper index NL denotes the nonlocal interaction vertex. It should be noted that to keep 
the U{1) 

cm gsuge invariance, a Wilson's line connecting the charged particles at different positions 
should be introduced. The kinematic parameter is defined as 

with Mi as the mass of i—th constituent. And correlation function ^(y^) describes the distribution 
of the constituents in the molecule. The Fourier transform of the correlation function reads 

Hy') = 1 0^Hp')e-"^-' (12) 

In the following numerical calculations, an explicit form of $(p^) is necessary. Throughout this 
paper, we will take the Gaussian form 

^p') = exp(pVA|) (13) 

where the parameter As is a size parameter which parametrizes the distribution of the constituents 
inside the molecule. The magnitude of is determined by requiring that the effective coupling 
constant should be stable against it. It should be noted that the choice (jl3p is not unique. In 
principle any choice of ^(p^), as long as it renders the integral convergent sufficiently fast in the 
ultraviolet region, is reasonable. In this sense, can be regarded as a regulator which makes 

the ultraviolet divergent integral well defined. In addition, we would like to point out that, in the 
limit ^{y'^) — > <J^(y), the interaction between the scalar meson and its constituents becomes a local 
one, i.e., Eq. ([TO]) approaches to Eq. dSj). 
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B. The Effective Coupling Constant Between the Scalar Meson and its Constituents 

After the discussion on the interaction Lagrangian of the scalar meson and its constituents, 
we are ready to calculate the effective coupling constant with the help of the compositeness 
condition 



Zs = 1- gl^'simh Ml M|) = 1 - gl^^siq^; Ml M|) 



(14) 



where Mi and M2 are the masses of constituents Pi and respectively, g^^s is the mass operator 



of scalar meson which is depicted in Fig. [2l This compositeness condition means, after renormal- 
ization, the scalar meson degree of freedoms are removed from the original (bare) Lagrangian and 
their dynamics are substituted by the relevant constituents. In addition, it also indicates, in this 

model, the scalar meson cannot arise as a final or initial state since, according to the LSZ reduction 

1/2 

rule, each such external state contributes a factor to the physical matrix element so this kind 
of matrix elements vanish. 




FIG. 2: Mass operator of scalar meson S. 

Using the compositeness condition (jl4p and the constituents contents of the scalar mesons, one 
can get the expressions of the coupling constants as 



1 

"0 



1 
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S„+(m„o;m^+,m^o) 

■'- ^ I 2 2 2\|V^//2 2 2 m 
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/ / 2 2 2\ 

^.+ (m^*;m^+,m^) 



g2 

^0 



sr-i t 2 2 2\ 
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Sin ^5 2C0S^ 2W ,2 .2 2s 



1 sine's , cos6ls.2rw/ 2 2 2 x , ^1 r 2 2 2 ^ 
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+ C^+'z^)^^U-nl;,n%n4) (15) 

The mass operator S'^(m|.; Mf, M|) can be calculated by evaluating the standard loop integral. 
From the diagram shown in Fig. [21 in the local interaction case we have 

TP,. 9 9 9s 1 rM?-M|, Ml 

S^i^ '(ml; M?, M^) = ^-^<^ ^ ^ ^ n — - 1 



mi.\/—X ^-^ a/— A J 



m^V — A 



where 



z± = m| ± (Mf - M|) 

A = A(m|, m2,m|) (17) 

with A as the Kiillen function. And similarly, for the nonlocal interaction case, the standard 
evaluation leads to 

1 l-CO l-OO p J 

SgL'(m|;Ml^M|) = ^ / dadB [-f-^\z^)] (18) 

16vr^ Jo Jo (1 + a + f3Y dzm 



where 



Ps = uJla + ujII3 + aP 



m| - aMl - pM^ (19) 



1 + a + P 

It should be noted that the expression (llSp is independent of the explicit form of <I>. And in the 
derivation of (jlSp we have applied the Laplace Transform 



Fis) 



/ f{s)e~'^dx] for Res > (20) 
Jo 



Up to now, the numerical results of the coupling constants , 5^0 , gj^*+ and g^*o can be yielded 
by inputting the physical masses of the relevant mesons. But the coupling constants gfg and g^ can 
not be got without determining the mixing angle. In the following, we will not discuss the coupling 
constant ga because of the large uncertainty of the sigma meson mass which should be used in 
the compositeness condition to determine the coupling constant g^. For other scalar mesons, the 
explicitly input masses are [l| (we adopt the central values of the scalar meson masses) 

map = 985.1 MeV; Mk* = 672 MeV; = 980 MeV 

mK+ = 493.677 MeV; m^o = 497.648 MeV; = 547.51 MeV (21) 
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With these numerical values, we get the following numerical results for the local interaction coupling 
constants 



^LC ^ 3 554 LC ^ 3 -^gg lc ^ ^^2.99 GeV ; g^^^ = 13.06 GeV 



9 



1 , , sin Oq cos 6*5 , 9 „ , cos Oq 2 sin 0o , o 
= 0.1194 X ( ^ + — + 0.0117 X ( — ^ + 



y/6 



V3 



V6 



(22) 



In the following, we will fix the coupling constant gf^ and mixing angle 6s using the two-photon 
decay of /n. The two-photon decay of /o 27 has been studied in the hadronic molecular model in 
Refs. 17,Q)[3l- the local interaction case, one should consider the diagrams shown in Fig. [3l 
And, even in the nonlocal interaction case, only including these diagrams is enough since, to our 
experiences, the contributions from diagrams with photon emerges from the nonlocal vertex are 
negligible. 





FIG. 3: Diagrams contribute to the decays of /o 27. 



Generally, the width for decay fo — > 27 can be expressed as 



r(5^ 27) = 167Tai^r,^G\mi)m's; 



(23) 



And, because there are only charged kaon mesons in the loop, g^^ relates to the coupling constant 
g^.^ via relations 



cos d sm U 

For the local interaction case, after standard calculation we get 

»1 rl-x^ 



1 



X1X2 



-xiX2m'j^ + m\ 



And for the nonlocal interaction, we get 

1 1 



167r2 A| 7o 



111 

da^da2da■^ 5 ( — \- ai)( — h 02)^(2/?) 



with 



zr 



11 1 1 ^ 

——3 (- + a2 + a3)(^ + a2)m\ + (- + a2)m% - ^ 



(24) 



(25) 



(26) 



(27) 
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Using the average data quoted by PDG ^ 

r(/o ^ 27) = 0.29+°:°^ KeV (28) 
and the expression (1250 . we have the numerical results of the coupling constant gf^ as 
4c = 2.878 GeV; 

g^^ = 2.941 - 2.915 GeV; for A5 = 3.0 - 4.0 GeV (29) 

which leads to 

g^^ = 4.175 GeV; for ^5 = 42° 

g^^ = 4.182 GeV; for ^5 = 68° (30) 

In principle, with the help of the decay of o" ^ 27 one can determine the exact value of Og if this 
model can explain all the data well. But, as we mentioned above, we will not study the physics 
of fj meson because of its large mass uncertainty, so that at present we will adopt 63 = 68°. We 
would like to point out that our final results do not depend on the choice of 9s closely, since the two 
angles lead to the same effective foKK coupling constant and the rjr] component plays negligible 
roles in the quantities we calculated. 

Using the same method, in the nonlocal interaction case, one can also determine the effective 
coupling constants g^. In this case, since the coupling constants are functions of A5, one should 
determine the magnitude of A5 at first. Our principle is that the effective coupling constants should 
be stable against Ag. With this criterion in mind, and running Ag from 1.0 GeV to 10.0 GeV, 
we find the physical region of A5 should be 3.0 — 4.0 GeV. We list the numerical results of g^ in 
Tablejl] and compare them with that from the local interaction case and other literature. 

In this table we find in the physical region of A5 all the coupling constants calculated in the 
nonlocal interaction case agree with that yielded in the local interaction case. However, compared 
with Ref. 1^, one find that our result for g^^ is much larger. This is because, in the present 
construction, in addition to the KK constituents, /o also consists of r]r] component so there is 
the mixing angle 63- When the mixing angle 63 is included our effective coupling constant Qf^^K 
agrees with that given in the literature. This is one of the typical properties of the present work. 
Concerning the typical property that numerical values of the coupling constants yielded from the 
nonlocal interaction case agree with that from the local interaction case, in the following calculation, 
we will adopt the local interaction vertex. 
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TABLE I: Numerical results of the coupling constants (in unit of GeV). The range of our results is due 



to the variation of A5 from 3.0 — 4.0 GeV. 







9< 


9k^^ 


9x^0 


5/0 


0s 


Our results (NL) 


3.591 - 3.577 


3.252 - 3.227 


14.13-13.73 


13.23- 13.81 


4.230 - 4.210 
4.237-4.217 


42° 
68° 


Our resuhs (LC) 


3.554 


3.169 


12.99 


13.06 


4.175 
4.182 


42° 
68° 


Ref. \18] 










3.09 




Ref. [19] 


3.61 (NL) 
3.45 (LC) 


3.61 (NL) 
3.45 (LC) 






3.06 (NL) 
2.87 (LC) 




Ref. [33] 










3.27 ±0.99 





At last we would like to mention that, because of the KK components in /o and aS, there is 



Og — /o mixing in the hadronic molecular model. In the present work we will not study the mixing 



effect of these two mesons. For persons who interest in this topic please see Refs. 
the relevant references therein. 



19|, 



35( 1 and 



C. The Leptonic Decays of the Scalar Mesons 

Next, we will calculate the scalar meson leptonic decay widths in the molecular picture. It 
is well known that, for a pseudoscalar meson, for example vr"'", the leptonic decay constant /^r is 
defined by 

(|M7^75d|7r"(p^)) = -ifnP^ (31) 

And since leptons are free from strong interaction, one can express the width of the pion weak decay 
7r~ lT>i in terms of the leptonic decay constant f-,^ using naive factorization. Comparing with the 
data, one can extract the numerical value of f-,^. On the other hand, if the wave function of pion is 
known, we can also calculate this quantity by the standard loop integral. For our present problem, 
the scalar mesons are interpreted as hadronic molecules and the coupling constant between the 
scalar meson and its constituents is determined by the compositeness condition. So that we can 
calculate the leptonic decay constants of scalar mesons via the standard loop integral. At the 
meson level, we define the leptonic decay constant of the scalar meson via 



{S{p)\Pid^P2\()) = ip^fs 



(32) 
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where is the momentum of scalar meson and Pi and P2 are the two constituents of scalar meson 
S. To calculated this quantity, one should concern the the Feynman diagram depicted in Fig. HI 




FIG. 4: Diagram relates to the leptonic decay constant of the charged scalar meson. 



The coupling constants between the weak gauge bosons and pseudoscalar mesons can be yielded 
by gauging the nonlinear sigma model and relating the flavor symmetry of the nonlinear sigma 



36| 



model to the flavor symmetry of QCD. The relevant terms are 

^gauged ^ ^Tr[Z)^[/(D^C/t) + xt[/ + C/tx] 



(33) 



where x = 2BM with i? as a constant relating to the quark condensation and M = 
diag(mu, mdjiris) as the quark mass matrix. The field U{x) is expressed in terms of pseudoscalar 
fields as 



Uix) = exp{i(f){x)/F^) 



/ vfO + ^r/ V2TT+ V2K+ \ 



(34) 



V2tt~ -^o + -i=r/ V2K^ 

where F^^ is the pion leptonic decay constant with the value F-,^ = 92.5 MeV. Under SU{3)l x 
SU{3)r chiral transformation meson matrix U{x) transforms as 



U{x) ^ gRUix)gl 



The covariant derivative is defined as 



Df,U = d^U-iAR,^U + tUA 



(35) 



(36) 



Al-^ and Ar-^ are the gauge fields corresponding to the gauged left- and right-handed chiral 
symmetry, respectively. 

By matching the chiral symmetry of QCD and the transformation of the field U{x) one can 
express these gauge fields in terms of the electroweak gauge bosons as js?! 

AR-y, = -eQA^ - g — QZ^ 

cos 6w 
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^L;m = -eQA^ + gQzZ^ + ■j={W;Qw + W-Q\^) 



(37) 



where Q is the charge matrix of quarks and in the three flavor case Q = diag(2/3, —1/3, —1/3) and 
e = gsinOw- The matrices Qw and Qz are defined as 

^ Vud Vus ^ 

Qz 



Q 



w 



us 






1 



cos 6 



1 1/2 





\ 





-1/2 





V 





-1/2 / 



sin 9w 



cos 9w 



Q 



(38) 



where Vud and V^s are the appropriate Cabibbo-Kobayashi-Maskawa matrix elements. And g is 
the coupling constant of the SU (2) l weak gauge group in the standard model and at the lowest 
order perturbation theory it is determined by the Fermi constant and the W boson mass via the 
relation [l| 

^2 



Gf = V2 
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1.16637(1) X lO^^GeV" 



Explicitly, one can get the following WPP interaction Lagrangian 



rW 



4^/2 



(39) 



(40) 



As an example, let us consider the decay of Cq (980) — > Ivi. Its matrix element, on the one 

hand, can be expressed in terms of the leptonic decay constant f + , and on the other hand, can be 

"o 

calculated directly from the relevant Feynman diagram in terms of the coupling constant o + , so 

Op 

that f + can be expressed as a function of the coupling constant o + . In this sense the numerical 
value of f + is calculable in the molecular model. 
Consider the following effective Hamiltonion 



= V2iGFVUK%K-][h^PL 



(41) 



on the one hand, the relevant matrix element, with the help of "naive factorization" , can be written 
as 



iM{a-^^lvi) = -V2iGFVud{hi>i\[K''d^K-][h^PLyi]\a^ 
= V2GFVudP^ f-u{p^)-f^PLv{pu^) 



(42) 



On the other hand, this matrix element can be expressed in terms of the coupling constant g^-\ 
and standard loop integral as 



iM{aQ Ivi) = -V^GFVuduiPf.h^PLvipu^) 
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/d^k I 1 

(27r)4 - m\, {k + vY - ^ ^'^o ^^^^ 

Comparing (I42p and (I43p and after some calculations, one gets the leptonic decay constant of ag 
as 

= 5„+-^("ia,p"i^o,m^+) (44) 

where 

/«,<o,<.) = /^x(l-2x)ln[ (^-^)"^^°+^"f -^(^-^^< j (45) 

where \x is the dimensional parameter introduced in dimensional regularization. It should be noted 
that, in this expression, the term (1 — 2x) In /i^ vanishes after the Feynman parameter integral. In 
this sense the the expression ()45p is scale independent. With the Eqs. (j44l45p we get the numerical 
result of the leptonic decay constant /^+ for scalar meson as 

/a+ = 5,+ ^L«,,m2,+ ,m^o) = 0.1530 MeV (46) 

Similarly, one can get the leptonic decay constant for K^^ as 

= 9kI+ {ml, ,ml+,ml) = 3.463 MeV (47) 

For Kq^ meson, one can also define its leptonic decay constant as (p2|) and express /^.o as (jH]). 
The numerical calculation yields 

Iks^ = dKS^himl, , m|,o, m^) = 3.208 MeV (48) 

From this result one can see that the isospin violating effect is important for the study of Kq 
mesons in the molecular model. 

We would like to mention that, in contrast to the quark model where the leptonic decay constants 
of neutral mesons are identified with their charged partners, the leptonic decay constants of neutral 
molecular scalar mesons /o and are zero since, equation (|45p vanishes in case of Mi = M2. 
Physically, this is because the weak interaction (and the electromagnetic interaction) is mediated 
by vector current of pseudoscalar mesons. 

With the help of leptonic decay constants calculated above, we can calculate the leptonic decay 
of the charged scalar mesons. In terms of the leptonic decay constants, one can express the partial 
width for decay Ivi as 



/ maomAl Y 

87r a+ \ mi 



r(ao- - m = -^^f\rna,mf[l--f (49) 

''ao ' 
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Along the same derivation, the leptonic decay width of Kq — > Ivi can be expressed as 

rW-'-,) = 5%^4»..-»?('-;|:)'^ (50) 

The numerical results are found to be 

T{a- e-pe) = 3.091 x 10"^° KeV ; T{a~ ii~v^) = 1.306 x 10~^^ KeV 
T{K*- e'Ve) = 1.620 x 10"^^ KeV ; r(K*- ^ ii~u^) = 6.755 x lO^^^ j^^y ^^^^ 

From these numerical results we conclude that, in the hadronic molecular interpretation, it is 
difficult to search for the leptonic decays of the charged scalar mesons in the near future. Or 
inversely, if the observed leptonic decay widths of scalar mesons are much larger than the present 
results, the hadronic molecular interpretation is suspectable. 



III. STRONG DECAYS OF LIGHT SCALAR MESONS 



In this section, based on the hadronic molecular explanation, we will study the strong decays of 
the light scalar mesons with masses below 1.0 GeV. Explicitly, we will study the scalar meson to 
two-pseudoscalar meson decays, i.e., the decays of —>■ rjTr~^,aQ —>■ 7]7r^,KQ~^ — > Ktt,Kq^ —>■ Kit 
and /o TTvr. These processes are important since they are the dominant channels for the relevant 

n 

scalar mesons [1]. To study these decays, one should consider the two kinds of Feynman diagrams 
depicted in Fig. O Diagram (A) is from the four-pseudoscalar meson vertex, while diagram (B) 
arises from the final state interaction. 



P2 




FIG. 5: Feynman diagrams contributing to the strong decay of ^ P1P2 {Mi denotes the constituent 
meson, Pi denotes the final pseudoscalar meson and V denotes the relevant vector meson). 



In the explicit calculation, we need the four-pseudoscalar interaction vertices which give diagram 
(A) of Fig. [5j Here we only consider the terms from the leading order chiral perturbation theory 



with the explicit chiral symmetry breaking terms 

i72 
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24F2 



TT 



+ 
+ 

+ 

+ 



+ 



(52) 



where only terms will be used in our calculation were written down. And to derive these relations 
we have used the Gell-Mann-Okubo mass relation 



4m^ = 3m^ + 



(53) 



The interaction Lagrangian for vector-pseudoscalar-pseudoscalar meson vertex can be written as 



(54) 



where the summation over isospin indices is understood and Ad^j^B = Adf_iB — Bdfj_A. The coupling 
constant gK*Kn can be fixed from the data for the strong decay width K* Kn. In particular 
the strong decay width T{K* — > Ktt) relates to gK*Kn via 

o2 



(55) 



where Pt^k* is the three- momentum of vr in the rest frame of K*. Using the data for the decay 
width one can deduce gK*K-K = 4.61 jl]. The coupling constant gx'Kq can be related to the gK'Kir 
using the unitary symmetry relation 



gK*Kri 



Fr, 



-gK'Kn = 6.14 



(56) 



Generally, the matrix elements corresponding to these two kinds of diagrams can be written as 
iM^'^\S ^ P1P2) = igsg,,^,,^p,pJ^^\ms;Mi,M2;mi,m2) 

iM^^\S^PiP2) = ig,gy,,^p^g^,,^^^I^^\ms;mv;M,,M2;mum2) (57) 
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where the notation rrii is the mass of the final state, 9]si-^m2PiP2 couphng constant of four- 



pseudoscalar meson vertex which was taken from the leading order Chiral perturbation theory 



39|] and 

9vM p is coupling constant of vector-pseudoscalar-pseudoscalar meson vertex. The 
functions I*^"^^ and I^^^ can be calculated using the standard technics of loop integral. Concerning 
the effective Lagrangian ([52]) . we will choose 5aj^m2PiP2 ~ ^/(^^-^w)- 
So we formally write the matrix element of 5 ^ P1P2 decay as 

iM{S^PiP2) = iM^^\S ^ P1P2) + iM^^\S ^ P1P2) (58) 

The decay width, in terms of I^^^ and can be expressed as 



r{s^PiP2) 



9 



^ ^ 2 2 2^ 



M1M2P1P2 



/(^)(m5;Mi,M2;mi,m2) 



+9vMi Pi 9vM2P2 ^^^^ ('^S ; Ml , M2 ; mv ; mi , m2) (59) 

where Pcm(m,|; mf, 771,2) ~ 2ms '^("^l' "^i' "^2) ^^^^ ^ Kallen function. 

Now, we will calculate the width for strong decay — > tt^t]. At first, let's consider the 
function J^"^). From the interaction Lagrangian given in Eq. (j52p . one can see that the four- 
pseudoscalar-meson vertex consists of two parts: the part includes derivative and the other part 
without derivative. So that we can formally do the decomposition 

l'-^\mao;rnK+,mKo;mrj,m^+) = 9p^P2PiPj''D\'^ao'^^K+,rniio-,mr^,m^+) 

+9p3'2P3Pa^n'd ("iao ; iT^K+ , m^o ; m^, ) (60) 

where the upper indices D and ND denote the contributions from Lagrangian with derivative and 
without derivative terms, respectively. After explicit calculation, we get 

4^^(m.o;m;,.,m^o;m„m..) = J (^•^(^') ^2 _'^2^ (k + - mj,, ^^'^ 

X [{k + p)+ pi] ■ (k - P2) + {pi - k) ■ [-pi -{k+ p)] 
J A) , ^ /" A^k i i 



I^j,{m.,;mj,,,m^o;mr„m^,) = J ^^(fc ) ^2 _ ^2^ + ^)2 _ ^2^^ (62) 

where p is the momentum of the incoming meson, and pi and P2 are momentum of the outgoing 
rj and vr"'", respectively. Since the momentum integral is divergent, the form factor J-{q^) was 
introduced to suppress the divergence. Explicitly, the momentum-dependence of the form factor is 

ne, ^ (63, 
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where n = 1,2 correspond to the monopole and dipole forms, respectively |40l |. Through out 
the following calculation we will select the dipole form since the above integral is quadratically 
divergent, i.e., n = 2. For the parameter A, since we only include the diagrams with the exchanged 
mass up to , we will take the typical value A = 1.0 GeV. In Appendix. [XI we express ll^^ and 
p in terms of the standarc 
For diagram (B) we have 



(A) 

Ij^l) in terms of the standard n-point scalar and tensor integrals 



l''^>imao;mK+,mKo;mr^,m^+) = / 77r^-^(^^)7] ^ 2~77~ vi 2—7:2 2" 



X (2pi - k)^ [g^, - A^] {k + 2p2)u (64) 

which is also expressed in terms of the standard n-point scalar and tensor integrals in Appendix. lAl 
Substituting the relevant masses and coupling constants and with the help of the software 
package Loop Tools 



4ll |. we get the width for the strong decay Oq TT~^r] as 

r(aj[ ^ vr+r?) = 58.81 MeV (65) 

Along the same method, one can get the following decay width 

r(a[] ^ 7r°r/) = 59.21 MeV (66) 

Similarly, for other strong decay widths we have 

r{K*+ K°TT+) = 11.06 MeV (67) 
T{Kf K-TT+) = 12.37 MeV (68) 
r(/o ^ TT^TT') = 30.65 MeV (69) 

It should be noted that the analytic forms of diagram (A) in Fig. [5] for these three processes are 
different from that of oq Tjn decay we listed in Appendix. [Al For simplicity, we will not write it 
down explicitly. 

With the help of the isospin relation, we get the following decay widths 
r(ir*+ ^ K-k) = r(K*+ ^ K^-n^) + r(i^*+ ^ i^°7r+) = ^r(i^"*+ K^tt-^) = I6.59 MeV 
r(i^o° ^ Ktt) = r(Ko° ^ K+TT') + r(i^o° ^ K^t^^) = ^r(i^o° ^ K+n') = 18.56 MeV 
r(/o ^ vr^) = r(/o ^ 7T+7T-) + ^r(/o ^ 7r°7r°) = 45.98 MeV (70) 

We summarize our results and compare them with that given in the literature in Table . HIl 
From this table we see our results for r(/o — > vrvr) and r(ao — > ryvr) agree with that of Ref. [191], 
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TABLE II: Numerical results of the widths for the strong decays S PP (in unit of MeV). 





Our results 


[1^ 


[8] 


[1^ 


[42] 


PDG [1] 


Meson structure 


Molecule 


Molecule 


hadronic 


qq 








45.98 


57.4 


18.2 


53 


58 - 136 


40 - 100 


a+ ^ n+ri 


58.81 


61.0 


21.0 


138 


98 


50 - 100 


a[j vrOr/ 


59.21 


61.0 


21.0 


138 


98 


50 - 100 


K*+ Kit 


16.59 




500.0 


193 


241 - 251 


550 ± 34 


Kn 


18.56 




500.0 


193 


241 - 251 


550 ± 34 



and the tiny differences can be understood by concerning tliat the scalar-pseudoscalar-pseudoscalar 
vertices used in the present work are local one but that in Ref. [3] are nonlocal. This indirectly 
indicates that the scale A = 1.0 GeV we chose is reasonable. For the /o vrvr decay, our result 
is consistent with that from both the qq and tetraquark interpretations and all the results are 
consistent with that given in PDG. For the oq rjir decays, the results from both the hadronic 
molecular interpretation and tetraquark interpretation are consistent with the data but hadronic 
interpretation leads to a smaller result than the tetraquark interpretation. At last, let us turn to 
the results of Kq decays. One can see that our result is much smaller than the data and other 
theoretical approaches. In fact, concerning our yielded result is based on the central value of 
Kq mass, rrix* = 672 MeV [ll], we varied rrix* up to 1.0 GeV as a check and it is found that 
T(Kq — > Kit) < 70 MeV which is still much smaller than the data and other approaches. In this 
sense, it is difficult to arrange the scalar mesons with masses below 1.0 GeV into the same nonet. 

At last, we would like to mention that, in contrast to the hadronic interpretation, the the large 
decay widths of Kq in tetraquark picture can be easily understood at the quark level. At the 
quark level, in the tetraquark interpretation Kq ~ [gs] (with q as unflavored quark), it can 
easily decay into Kir by interchanging a pair of quark and antiquark and this process is OZI rule 
superallowed. However, in the hadronic molecular model, because of the ss component of rj meson, 
Kq ~ Kr] ~ [gs][ss], this process happens via an annihilation of a strange quark and antistrange 
quark, with subsequent qq creation, i.e., [9s][ss] — > [qgs] [qq] + [qs], so this process is OZI rule 
sub dominant. 
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IV. THE PRODUCTIONS OF NEUTRAL SCALAR MESONS /„ AND ag IN THE 

RADIATIVE DECAYS OF </> MESON 



In this section, we will study the productions of the scalar mesons /o and oq in the radiative 
decay of vector meson (p with including the intermediate axial-vector mesons and rji] component 
of /o meson. These processes are important because they have long been accepted as a potential 
route to reveal the natures of scalar mesons /o and ag. To study these processes, three classes 
of diagrams shown in Figs. [U [7] and [8] should be included. In Fig. [6l diagram (A) arises from 
the gauge of the derivative coupling of vector-pseudoscalar-pseudoscalar meson, while diagrams 
(B) and (C) are from the gauge of kinetic terms of the charged constituents. The diagrams in 
Fig. [7] contribute to both /o and oq productions from the intermediate axial-vector mesons and 
the Fig. [8] only contributes to /o production due to its 7]r] component. It should be noted that the 



diagrams like that in Fig. [7] but substituting the axial- vector mesons with vector mesons wi 
be considered in the present work since the sum of the (p and uj contributions almost cancel 



1 not 



43|. 




(A) 



fo/a« 




/()/«() 





/ll/ 0(1 



FIG. 6: Diagrams contributing to the (j) — s- /o/ao7 Decay (without intermediate resonance contribution). 




FIG. 7: VMD diagrams contributing to the cj) ^ /o/ao7 Decay. 

Before explicit calculation, we will discuss the interaction Lagrangian to be used in the following. 
At first, the (pKK interaction Lagrangian can be effectively written as 

^4>KK = ig4,KK(pfiK\x)d^K{x) (71) 

Using the data for the strong decay width r(</) K+K' = 4.26 MeV x 49.2% = 2.096 MeV [l|, 
one can fix the coupling constant g^pKR- In particular the strong decay width T{(j) — > K^K~) 
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FIG. 8: VMD diagrams contributing to the </) — > /07 Decay (contribution from the 77 constituent of /o). 



relates to q^kk via 



rf 



94>KK d3 



Gvrm 



— P 



(72) 



where P,^/^ is the three-momentum of K in the rest frame of (p. Using the masses for the relevant 
mesons one can deduce g^xK = 4.48. 



The interaction Lagrangian between photon and vector meson is given by 

. .,-2. 

?2. 



4i 



■■V-y 



-4F^egp^^Af,Ti[QV^ 



-A 



(73) 



3 ^ ^ 3 

where Q = (2/3, —1/3, —1/3) is the quark electric charge matrix and, the vector meson matrix is 
chosen as 

= \ V2p- -Z + c^ V2K*^ (74) 

and , ^. = 5.98. the u„.™™l,.. coupling—. 

To include the axial- vector mesons, according to PDG [1[, one should consider two nonets with 
quantum numbers J^*-" = and 1++. Explicitly, they can be written in U{?>) matrix forms as 

/a? + /i(1285) V2at V2K+j^\ 



A (1285) y/2Klj^ 



(75) 



\/2a^ —c 
\ V2K-^ V2K% /i(1420) J 

^6? + /ii(1170) V2bl V2K^j^ ^ 

V2b^ -6? + /ii(1170) V2K^B 
\ V2K-^ ^i(1380) j ^ 

where A^ and are the axial- vector matrices with quantum numbers J^*-" = 1^+ and 1"' , 
respectively. The mixture of Kia and Kib with approximate 45° mixing angle gives the physical 
states Ki{mQ) and Ki{Um). Explicitly 

i^i(1270) = -^{KiB-iKiA) 
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i^i(1400) = -^{K^B+iKiA) 



(76) 



Due to the C parity, the interaction Lagrangian for the axial-vector-vector-pseudoscalar cou- 
phngs have the following forms 



^AVP = igAvp^ 



^BVP = ^BVpTr B^{V^,P} 



(77) 



where i in front of g^^yp was added to keep the Hermitian of the Lagrangian. Using the concrete 
matrix forms, one can write down the interaction vertices that we are interested in 



+V2ff, 



K+^(1400)i^- - i^+^(1270)K- + KO^(1400)i^o - Ki%(1270)i?o 

Cbvp = V2g^yp [i^+^(1400)p°K- + K+jmO)plK- 

+V2K+^{im)p-K'' + ^/2K+^(1270)p;^o' 



+ H.C. 



-V2g 



K+^{im)K^ + K+^{1270)K- + i^?.^(1400)i?° + Kf.^{1270)K^ 



i^<AM^i;M(1380)r? + ^^c^;.(6° ^ + /ii;;.(1170)ry 



+ H.C. 
(78) 



The coupling constant g^^ can be determined by the decays of Ki — > pK via the expression 



3 H T^lciV 



(79) 



where g^^yp is the relevant coupling constant and Px^y is the three momentum of meson in 



the rest frame of Xi— meson. Explicitly, g 



K]^(1270) + p0_fs-- 



V2{g 

AVP + 9bvp) and g 



^^{9 EVP ~ 9avp)- Using the central values of the branching ratio and total widths from PDG 
r(A"i(1270) ^ pK) = 37.8 MeV and r(A"i(1400) ^ pK) = 5.22 MeV one can get 
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K-i(127Q) + p^K- 



3.52 GeV; g 



Xl{1400) + pOA'- 



0.54 GeV 



(80) 



which lead to 



^Avp 



1.05 GeV; g. 



1.44 GeV 



(81) 
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We would like to mention that, in the following calculation, we will not include the diagrams 
arise from the (p — uj mixing. This is because, the interaction Lagrangian describing the (p — u) 



mixing IS 



451 



Q4>^.^^. (82) 



where the coupling constant Q is determined by the relation 

e 



46| 



- ml 



0.059 ± 0.004 (83) 



all the diagrams with (f) — uj are suppressed by factor e compared with the diagrams without (p — oj 
mixing. 

With the above discussions, we can calculate the decay width now. Generally, due to the gauge 
invariance, one can write the total matrix element as 

iM{(l){p) ^ S{p )j{q)) = iei,{p)el{q){gfj,up ■ q - Qf^Pu) e G^sj (84) 
and the effective coupling constant Gtj,s-y consists of two parts 

Gcf>s-r = ^057 + G'0S'7 (85) 
where C^^^ is from Fig. [6l and C^^^ is from Figs. [7] and [HI The decay width can be expressed as 



r(0-57) = ^\G^s,\'P;' (86) 



where P* = (m| — vn?g) / (2'm^) is the three-momentum of the decay products. 

With these discussions and selecting A = 1.5 GeV which means the resonances with masses 
below 1.5 GeV were included, one can do the numerical calculation. Using the effective coupling 
constants, we get the decay widths from the contact diagrams and the resonances as 

r((^ ^ klY = 3.080 X 10"^ MeV ; ao7)^ = 2.329 x 10"^ MeV 

r((A ^ klY = 6.670 X 10"^^ MeV ; r(0 ^ 007)'' = 3.353 x 10"^^ MeV (87) 

From these results we see, compared with the contact diagrams, the contribution from the axial- 
vector resonances is negligible. In summary we have the total decay widths and the branching 
ratios (using the total width F,^ = 4.46 MeV [l|) as 

^{(t> hi) = 3.081 X 10"^ MeV; T{(t) 097) = 2.329 x 10"^ MeV 
Br (0 ^ /07) = 6.907 X 10"^ Br (c^ ^ 097) = 5.223 x 10"^ (88) 
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TABLE III: The branch ratio of the — * fo/aQ-f decays. 



Decay modes 


Present resuhs 


SND 


CLOE 


CMD-2 


Br(0 -> /07) 


6.91 X 10-^ 


3.5 X 10-'* [47] 


(4.47 ±0.21) X 10-4[49] 


2.90 X 10-'*[51] 


Br((/) 007) 


5.22 X 10-^ 


8.8 X 10-5 [48] 


(7.4 ±0.7) X 10-5 [50] 


9.0 X 10-5 [51] 


Br(0^/o7) 
Br(0^ao7) 


1.32 


3.98 


6.04 


3.22 



and the ratio 

R = lit ^ = 1-32 (89) 

r(0^ao7) ^ ^ 

In Table. Hill we compare our results with that from the data. From this table one can see that, 
for the decay (p clq'j, the present branching ratio is consistent with the observed data while, that 
for the decay (j) /07 is much smaller than data. This fact makes R = Br((/) /o7)/Br(0 
aQj) = 1.32 disagree with the data. Actually, as mentioned in Ref. [3;], this is one of the main 
problems in the hadronic interpretation of the light scalar mesons. Another observation is that, 
the final state interactio n g ives a negligible contribution to these decays which agrees with the 



conclusion given in Ref. 
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V. DISCUSSIONS AND CONCLUSIONS 



In this paper, in the framework of effective Lagrangian approach, we studied the properties of 
the light scalar mesons with masses below 1.0 GeV in the hadronic molecular interpretation. 

To determine the coupling constant between the scalar molecule and its constituents we 
applied the compositeness condition which has been used in our previous works. We found the 
numerical results of are consistent with that given in the literature. With the yielded coupling 
constant g^ we calculated the leptonic decay constants and leptonic decay widths of scalar mesons. 
Our numerical results show that the leptonic decays of scalar mesons are not observable in the near 
future experiments. 

The calculations of the strong decays conclude that the decay widths for /o tttt and oq r]Tr 
is consistent with the observed data while width for Kq — > Kir decay is much smaller than the 
data, even the ambiguity from the mass is considered. This observation shows that the hadronic 
molecular interpretation of Kq is disputable and, concerning the /o — > vrvr and oq — > tjtt decays are 
consistent with the data, the classification of the scalar mesons with masses below 1.0 GeV to form 
a nonet is disputable. 
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To study the productions of /o and oq in the radiative decays of </> meson, we included the 
contributions from the final state interaction, i.e., the contributions from the intermediate axial- 
vector mesons. The explicit calculation shows that the dominant contribution is from the contact 
diagrams of kaon loops and axial-vector mesons play a negligible role in these processes. The 
branching ratio for </> — > ao7 is consistent with the data while that for (j) /07 is much smaller 
than the data. This is another problem of the hadronic molecular model. 

We would like to say that, if the isosinglet pseudoscalar constituent r] is substituted by r]' , the 
width for decay Kq — > Ktt is not improved but suppressed. In fact, when rj meson is substituted 
by 7]' meson, the coupling constant g is improved to g ,± = 20.54 GeV, but coupling constant 



is suppressed to gK*Kri' = 1-095 when the mixing angle Bp = —9.95° 



52[ | and fg 



3K*K-q' 

1.26/7r [53| are applied. So that, compared to the rj meson case, the partial width T{K'^'^ Kir) 
from the final state interaction is suppressed by a factor (20.54 x 1.095)^/(12.99 x 6.14)^ ~ 0.078. 
For the contribution from the contact diagram, explicit calculation shows that the partial width 
due to this diagram is about 34% of that of the r] meson case. So that the partial width from the 
sum of these two diagrams is about one third of that of the r] meson case. 

One may notice that another possibility to improve the numerical value of the width for decay 
Kq Ktt is to enlarge the parameter A in the form factor !F{q^). In fact, we check that to make 
the numerical value of the width for decay Ktt agrees with the data, A ~ 1.5 GeV. But in 

this case, ^{clq — > rjir^) = 749 MeV and r(aQ — > rjir^) = 665 MeV. Both of them are much larger 
than the data. 

At last, we would like to mention that, we have estimated the decay of ^0(1430) Kir rudely. 
In this case, the singlet must be ij' . Our result for the partial width is ~ 124.0 MeV which is about 
a half of the data and the improvement is mainly from the phase space. In this sense, it seems that 
(T (600), /o (980), 0(980) and Kq{1A30) can be classified into a same nonet. However, this deserves 
further systematically study. 

In conclusion, our explicit result shows that, in the hadronic molecular model, it is difficult to 
arrange the scalar mesons with masses below 1.0 GeV in the same nonet. 
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APPENDIX A: INTEGRAL FORMULAS 

In this appendix, we will derive some integral formulae of the one-loop integrals for ao — ^ r/Tr 
decay. To derive the relations we have used p = Pi + P2 and p^ = m^jpf = mf. 

1. Integral formulae for the contact diagram 

In this subsection, we will reduce the momentum integral for the contact diagram to the standard 

(A) 

loop integral function. For , one has 



[(k+p) +pi] ■ {k-p2) + (pi - k) ■ [-pi - {k + p)] 



-WP1P2P3P4 



- (3m| - m? - mi - - M^) 

xIo{ms;Mc„M^^;A) 

+4'^^ {ms; , Me, ; A) + ^ {M,, ; A) (Al) 



where the relevant functions are defined as 

1 r d^k 1 P-A^ 



A;2 - Ml {k + pf - Ml 



1 

m 



(K - Y^[^o(m|, Ml, Ml) - Bo{ml A^, )] 



fe2 - Ml 



= j^lAoiMl) - Ao{A') - {Ml - A^)A_Ao{A')] 

- ^ :[Boiml,Ml,Ml)-BoimlAlM^.^' 



'dA2 



-{Ml-A')-^Bo{mlA',Ml)] (A2) 



And for /J^2 we have 



= -ig^^p^p^p^ Iq {ms ; Mc, , Mc, ; A) (A3) 
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2. Integral formulae for the final state interaction diagram 

In this subsection, we will reduce the momentum integral for the final state interaction diagram 
to the standard loop integral function. After some algebra, one has 

l(^Hms;M^„M^,;mv;m„m,) = / (fc _ ^^)2 _ + ^^)2 _ ^2^ p T^^^ 

k k 

x{2pi - k)f,[g^^ - -^]{k + 2p2)y 
= Io'\ms; Mc^, Mc2;mv;mi,m2; A) + ^4"^ (my; A) 



rriy 



1 



-[1- ^(M^, -mi)]4')(Me,;my;m2;A) 
niy 

- [1 - ^ - (M,, , my ; mi ; A) 

my 

- + M^^ - m^ - 2m| + mf + m^ 

-^(Ml-mj)iMl-ml) 
rriy J 

x/_2(m5;Mci,Mc2;my;mi,m2; A) (A4) 

where 

( ) 1 f d'^k 

lo (ms; Mci, Mc2;mv;mi,m2; A) = - t^-TaH^'^) 

1 1 e-A^ 

^{k-pif- Ml {k + p2)2 - Ml P - m^y 

= {my - A'^)j^[Co{ml,ml,ms,Ml,my,mlj 

-Co (ml , mi , m| , , A^ , m^ J] 
7o (M.,;my;m2;A) = - J ^^(fc ) + ^^)2 , fc2 _ ^2^ 



= Y^[-^o(?«i,my,m^J - 5o(m|, A2,m^J 
-(m^ - A2)^Bo(mi, A2,m^J] 

4^)(M.,;my;mi;A) = ]/ ^.Hk\ A^) ^ J. _ 



my 



! [-Bo (jTi? , my , m^.^ ) - Bo {m\ , A^ , m^^ ) 



167r2 

-(m^ - A2)^5o(m?, A^ m^^. 
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I-2ims; , Mc2;mv;mi,m2;A) 



1 



{k - pif - Ml {k + p2f - Ml fc2 



167r2 



[Co {m\ , ml , m| , Ml ,my,m] 



"02 J 



-Co{ml,ml,ms, Ml , m; 



"02 } 



-{ml 



dA2 



Co (ml ,mlml, Ml , , ml )(]A5) 
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